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Magnetic spin excitations in Mn doped GaAs: A model study 
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Abstract. - We provide a quantitative theoretical model study of the dynamical magnetic prop- 
erties of optimally annealed Gai-^Mn^As. This model has already been shown to reproduce 
accurately the Curie temperatures for Gai-^Mn^As. Here we show that the calculated spin stiff- 
ness are in excellent agreement with those which were obtained from ab-initio based studies. In 
addition, an overall good agreement is also found with available experimental data. We have also 
evaluated the magnon density of states and the typical density of states from which the "mobility 
edge", separating the extended from localized magnon states, was determined. The power of the 
model lies in its ability to be generalized for a broad class of diluted magnetic semiconductor 
materials, thus it bridges the gap between first principle calculations and model based studies. 
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The prospect of manipulating the electronic spin for 
spintronics applications has generated a tremendous in- 
terest in the so called diluted magnetic semiconductors 
(DMS) over the past few years [HE]- Among the III- 
V semiconductors, Mn doped GaAs is the most widely 
studied. The understanding of its fundamental physical 
properties involves large theoretical speculations [5H3 , like 
there is still a controversy regarding the existence of a 
preformed impurity band (IB) in GaMnAs. Most of the 
model studies are based on a perturbative treatment of 
Mn in GaAs (the valence band (VB) picture) 03]. This is 
in contradiction to the ab-initio calculations [SJE]j which 
shows that Mn strongly affects the nature of the states 
close to the Fermi level (Ep) leading at low impurity con- 
centration to a preformed impurity band. This supports 
the IB picture and rules out the VB picture. Note also that 
(Ga,Mn)As is close to the metal insulator transition. In- 
deed after annealing as grown samples an insulator-metal 
transition was also often observed experimentally [8 , 9 . 
This important feature cannot be captured in the frame- 
work of the perturbative VB theory. 

Infrared and optical spectroscopy measurements have 
also shown that Ep resides indeed in a Mn induced IB 
in GaMnAs [ID] . In the most frequently studied model, 
based on a six or eight band Kohn-Luttingcr Hamiltonian 
[2J[3], the pd-coupling is treated perturbatively and the 
dilution effects are neglected (the VB picture), which as 
aforementioned are inconsistent with the ab-initio based 



studies. It has been demonstrated that the perturbative 
treatment is inappropriate for DMS and that both ther- 
mal fluctuations and disorder effects play a crucial role. 
The importance of treating disorder effects in a reliable 
manner was often underlined [1 1H13J . The first principle 
based calculations were able to reproduce the Curie tem- 
peratures accurately [blHT5] but they are essentially ma- 
terial specific. Hence the ideal tool to identify and analyze 
the effect of the relevant physical parameters is a minimal 
model approach based on a non-perturbative treatment of 
the substitution effects. 

In this communication we use the one band V-J model 
[13] . treated non-perturbatively, to study the dynamical 
magnetic properties of the III-V DMS GaMnAs. It will 
be shown that the V-J model provides the missing link 
between ab-initio and model studies. To be more specific 
we adopt a two-step approach to calculate the magnetic 
properties of GaMnAs. First the couplings between the 
magnetic impurities are calculated within the one band 
V-J model without the use of any effective medium the- 
ory. In the second step the effective dilute Heisenbcrg 
Hamiltonian is treated within the self consistent local RPA 
(SC-LRPA) theory to calculate the magnetic properties as 
well as the T^s. Note that the accuracy and reliability 
of SC-LRPA to treat dilution/disorder and thermal fluc- 
tuations has been proved time and again [T|)][T7]. The 
non-perturbative treatment of the V-J model was shown 
to evaluate qualitatively the essential magnetic proper- 
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ties of a whole class of III-V materials [IS] . Recently the 
model was also shown to successfully reproduce and ex- 
plain quantitatively the measured optical conductivity in 
Gai-^Mn^As [TS]. In the present work, wc calculate the 
Tcs, the magnon density of states (DOS) (average and 
typical) along with the "mobility edge" , the magnon spec- 
tral function and the spin stiffness as a function of the im- 
purity concentration in optimally annealed Gai_. r Mn r As. 
The term optimally annealed implies that the concentra- 
tion of compensating defects (As anti-sites or Mn inter- 
stitials) is negligible. Note that here we focus on homo- 
geneously diluted systems with no correlations in disorder 
(absence of inhomogeneities). 
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Fig. 1: (Color online) Curie temperature (in K) as a function 
of the Mn concentration (x). The squares represent our model 
and the circles correspond to ab-initio calculations. Inset (a) 
shows the MF-VCA Tc (in K) as a function of x. Inset (b) 
shows the exchange couplings (in meV) as a function of dis- 
tance for three different values of x. 

The one band V-J Hamiltonian describing the interac- 
tion between the carriers (holes or electrons) and the lo- 
calized impurity spins is given by, 
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where the hopping term ty =t for i and j nearest neigh- 
bors, otherwise zero. c\ a (c; CT ) is the creation (annihila- 
tion) operator of a hole of spin a at site i. Ji is the p-d 
coupling (J p d) between localized Mn spin Si (|S;|=5/2) 
and a spin carrier Sj (p-band). The on-site potential V* 
results from the substitution of Ga 3+ by Mn 2+ . Note that 
V, is directly related to the position of the bound state 
with respect to the VB. J;=PiJ p( z and Vi=p^V, where Pi=l 
if the site is occupied by an impurity, otherwise zero. Here 
x and p represent the Mn concentration and hole density 
respectively. The next task is to fix the model param- 
eters. Since our calculations are performed on a simple 
cubic lattice (1 atom per unit cell) the value of the lattice 
parameter used here is a=^jj=3.55K, where ao=5.65Ais 
the lattice constant for zinc-blende GaAs (4 Ga per unit 



cell). The hole effective mass for density of states in GaAs 
is known to be ~ 0.5m e [15], from which the value of t is 
fixed to 0.7 eV. Note that a variation of about ±10% of t 
docs not affect our results considerably. J p <zS is set to the 
value of 3 eV since J p d ~1.2 eV [5D] is already known for 
GaMnAs. The last parameter, the on-site potential V is 
the crucial one. It is set to 1.8t in order to reproduce the 
bound hybridized pd-states energy Eh psIIO meV in GaAs 
host [5T] with respect to the top of the VB. Now in the 
realistic material each Mn 2+ provides one hole and n;=3 
degenerate p-d states near the top of the VB (here p=x for 
the well annealed case). Thus for reasons of consistency 
our model calculations for optimally annealed samples are 
performed for the hole density p=p/n/. Then in both cases 
the IB is one-third filled (for more details see Ref. (T51l2"2"] y 
Let us now proceed with the calculation of the exchange 
couplings from our one band model. For a given impu- 
rity concentration l x' and disorder configuration 'c', the 
Hamiltonian (I) is diagonalizcd exactly in both spin sec- 
tors. We have performed the calculations on simple cubic 
systems with sizes varying from L=I6 to L=24 and the 
average over disorder is done for a few hundred configu- 
rations. A careful study of the finite size effects on the 
magnetic couplings has been made. The diagonalization 
provides us with the eigenvalues and eigenvectors denoted 
by |^ Q )}, where <r=t,i and a=l,2,...,N (N=L 3 is 

the total number of sites) , which are used to evaluate the 
couplings. The magnetic coupling between two localized 
spins is given by the generalized susceptibility |23| 
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where the Green's functions are defined as Gfj(u)= 
(ia\ — A— \ja). It is found that the couplings calcu- 
lated for well annealed Gai-^Mn^As are rather short 
range and essentially ferromagnetic. These are similar 
in nature to those obtained from first principle studies 
[T31 [H]. Now the couplings that are used to calculate 
the Tcs and other magnetic properties are defined as 
Jij(x,p)=niJij(x,p = p/n{), since we stress again that 
in our one band model each Mn 2+ provides a single state 
unlike the 3 p-d states in the realistic case. The exchange 
couplings for three different Mn concentrations are shown 
as a function of distance in inset (b) of fig. [TJ Note that the 
same argument was used to calculate the transport prop- 
erties, namely the optical conductivity, in Gai_ x Mn x As 
and excellent agreement with experimental results was ob- 
tained [IS]. As we have found that the spin stiffness is 
sensitive to the size of the lattice on which the couplings 
are calculated, here the calculations have been performed 
systematically on considerably larger systems compared 
to previous studies |22| leading to an improvement in the 
couplings at large distances. It should be noted that calcu- 
lating the exchange couplings and the corresponding mag- 
netic quantities simultaneously for a given configuration 
and then averaging these obtained quantities over several 
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configurations leads to the same results as the ones shown 
here. 
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Fig. 2: (Color online) Magnon density of states (along the y 
axis) as a function of the energy cj (in meV) (along the x axis) 
for different Mn concentration x. The solid (dashed) line rep- 
resents the typical (average) density of states. 

In fig. [T]we show the Curie temperature (in Kelvin) as 
a function of the Mn concentration for well annealed sam- 
ples. To evaluate the Tc, the effective dilute Heisenberg 
Hamiltonian Hn e i S =— . piPjJij(x,p)Si ■ Sj, is treated 
within the SC-LRPA theory. After performing the Tyab- 
likov decoupling and diagonalizing the above Hamiltonian 
we obtain the retarded Green's functions given by 
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where <S> c a and w£ are now the eigenvalues and eigenvec- 
tors respectively of Hneis for a given configuration 'c', and 
(S") is the local magnetization (for details see Ref. |16j). 
The diagonalization of the above Hamiltonian was per- 
formed on very large systems (typically L=50) together 
with a systematic average over a few hundred configura- 
tions of disorder. We have found that the Tcs shown here 
are comparable to those shown in Ref. [22] ■ As is clearly 
seen from the figure, our model calculations are in very 
good agreement with the Tps obtained from first principle 
studies jT5ll6l25| , which in turn reproduced accurately the 
experimental data [2"5H2"8] . We again stress the fact that 
there arc no adjustable parameters in our model calcula- 
tions. It is worth noting that this model can also be used 
to evaluate the Tcs and the magnetic and transport prop- 
erties for a whole class of DMS materials, which makes it 
all the more powerful [22]. Several attempts have been 
made to calculate Tcs using mean field theories [29] but 
they have always resulted in ovcrcstimations. As shown 
the mean field Virtual Crystal Approximation (MF-VCA) 
Tc (inset (a) of fig. [1} already leads to room temperature 
for only 2% of Mn. Once again it shows the clear overesti- 



mation in the Tc values and hence the crucial importance 
of thermal fluctuations and dilution effects. 

The average and typical magnon DOS as a func- 
tion of the energy for different concentrations of Mn 
are plotted in fig. [5] The size of the simple cubic 
system in this case is 44. The average magnon 
DOS is given by p a vg(v)=((l/Ni mp )YiiPi(w))c, where 
Pi(cj)=—(l/2Tv(Sf))ImQii(uj) is the local magnon DOS, 
Ni mp is the total number of magnetic impurities and (...) c 
denotes the average over disorder configurations (for de- 
tails see Ref. [T7]). As seen from the figure, the shape of 
Pavg(w) exhibits a multipeak structure and changes sig- 
nificantly with increasing x. At sufficiently high dilu- 
tion, close to the percolation, there is a drastic increase 
in the low energy DOS. This is due to formation of iso- 
lated impurity clusters at low concentrations which have 
their own zero-energy eigenmodes and these contribute to 
the average magnon DOS. However p a vg(^>) is not able to 
predict the nature of the magnon modes, whether they 
are extended or localized. For this we calculate the typ- 
ical magnon DOS [30] which is given by p ty p(uj)=exp((\n 
Pi(u)) c ). The typical magnon DOS is a local quantity and 
unlike p avg (iS) it provides direct access to the "mobility 
edge" separating the localized modes from the extended 
ones. The solid (red) line in the figure corresponds to 
Pty P (uj) and apparently there is no significant change in 
the overall shape with variation in dilution. We have eval- 
uated from p t y P (uj) the energy E c , which separates the 
localized from the extended states. The values are indi- 
cated in the figure. As we increase the Mn concentration 
E c varies from 8 meV for x=0.02 to 32 meV for £=0.08. 
This shows that even for relatively larger concentrations 
most of the magnon excitations consists of localized states 
(fractons) [3"T] . 

In the next figure, fig. [3] , we show the magnon spectral 
function A(q, w)=-(l/7r((S 2 )))Im<7(q, lu), in the (q,w) 
plane over the entire Brillouin zone corresponding to three 
different concentrations of Mn. The spectral function is 
directly accessible by inelastic neutron scattering experi- 
ments and provides direct insight into the magnetic exci- 
tation spectrum (details can be found in Ref. Q7]). The 
average was done over a few hundred configurations of dis- 
order but we found that increasing the number of configu- 
rations beyond 50 left the results almost unaffected. Now 
similar results were obtained for the excitation spectrum 
of Gai-^Mn^As [32] but it is to be noted the couplings 
in that case were calculated from first principle Tight 
Binding Linear Muffin Tin Orbital approach. We observe 
that well defined excitations exist only in a restricted re- 
gion of the Brillouin zone essentially around the T point 
[q=(000)]. The spectrum shows a significant broadening 
as we move away from the T point. This is in agreement 
with the typical DOS calculations shown in the previous 
figure. Thus we can see that our results are in very good 
agreement with those obtained from first principle studies 
cited above. 

Now we proceed further and calculate the spin stiffness 
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Fig. 3: (Color online) Spectral function A(q, u) in the (q, ui) 
plane for different Mn concentration (x). The energy axis (y 
axis) is in meV. (The size of the simple cubic lattice is L=44). 



as a function of the concentration of Mn. As a first step 
the magnon energy w(q) is plotted as a function of q 2 for 
different shown in fig.|U Here w(q) is extracted from 
the first peak of A(q, tu) in the (100) direction. The sys- 
tem size varies from L=2>2 to L=44. The slope of these 
curves gives the spin stiffness D(x) for various x. Note 
that the perfect quadratic nature of the dispersion curves 
supports our claim that the average over disorder config- 
urations is sufficient. In fig. [5] we have plotted the spin 
stiffness D{x) (in meV.A 2 ) as a function of Mn density. 
As stated before the lattice parameter used in our simple 
cubic system is a=^f73=3.55A, which is also used to cal- 
culate D(x) here. It is worth noting that using the zinc 
blende lattice parameter (ao) here would lead to stiffness 
values ~2.5 times larger than the present ones. Let us 
now discuss the results. We have found that D(x) is sen- 
sitive to the couplings at large distances which appear to 
be strongly finite size dependent. Thus, the calculation 
of D(x) is performed using the couplings from different 
system sizes ranging from L=\Q to L=24, leading to error 
bars for D(x) shown by the shaded region in the figure. 
It is to be noted that these couplings do not have a sig- 
nificant effect on the Tps. In the figure we have shown 
the values of D{x) obtained from ab-initio studies [52] . 
Note that the stiffness values shown in Fig. 4 of Rcf. [3"2"] 
actually correspond to DS where S=5/2 [33]. Available 
experimental data for well annealed and as-grown sam- 
ples [341l35| are also shown in the figure. Surprisingly we 
find that the spin stiffness values from our model are in 
very good agreement, for the overall range of Mn concen- 
tration, with those obtained from first principle studies. 
We also note that the agreement with the experimental 
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Fig. 4: (Color online) Magnon energy w(q) (in meV) as a func- 
tion of q 2 (in A~ 2 ) for different concentration of Mn (x). (L is 
the size of the simple cubic lattice). 



results of Sperl et al. [33] is for the well annealed sam- 
ple of thickness 200nm. This is consistent with the fact 
that our couplings are calculated for optimally annealed 
samples. In Ref. [33] the authors have reported a value of 
DpsIOO meV.A 2 measured by ferromagnetic resonance in 
as-grown sample for x—0.05, which is also in good agree- 
ment with our results but then the error bar is about 40% 
and in all likelihood this stiffness should increase if the 
sample is optimally annealed. The deviation between our 
results and some experimental data could be explained by 
the fact that our calculations are performed for optimally 
annealed samples. 

Thus to conclude, we have provided a quantitative 
model study of the magnetic excitations in Mn doped 
GaAs. The model parameters were fixed from the first 
principle calculations. We have calculated the magnon 
DOS, the"mobility edge" separating the localized states 
from extended ones, the dynamical spectral function, the 
spin stiffness and the Tcs as a function of the Mn concen- 
tration. We found a remarkable agreement between our 
theory and first principle studies for both the spin stiff- 
ness as well as the T^s. At the same time we were able to 
reproduce most of the available experimental results. The 
strength of the model lies in the fact that it can be applied 
to a whole family of DMS materials and this proves to be 
the link between model approaches and first principle cal- 
culations. 
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